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Abstract 

A  class  of  non-oscillatory  numerical  methods  for  solving  nonlinear  scalar  con¬ 
servation  laws  in  one  space  dimension  is  considered.  This  class  of  methods  contains 
the  classical  Lax- Friedrichs  and  the  second  order  Nessyahu-Tadmor  scheme.  In  the 
case  of  linear  flux,  new  I2  stability  results  and  error  estimates  for  the  methods  are 
proved.  Numerical  experiments  confirm  that  these  methods  are  one-sided  I2  stable 
for  convex  flux  instead  of  the  usual  Lip+  stability. 


1  Introduction 

We  are  interested  in  the  scalar  hyperbolic  conservation  law 

m  (ut  +  f(u)x  =  0,  (x,t)  E  R  x  (0, 00), 

'  \  u(x,  0)  =  u°( x),  16I, 

where  /  is  a  given  flux  function.  In  recent  years,  there  has  been  enormous  activity  in  the 
development  of  the  mathematical  theory  and  in  the  construction  of  numerical  methods 
for  (1).  Even  though  the  existence- uniqueness  theory  of  weak  solutions  is  complete,  there 
are  many  numerically  efficient  methods  for  which  the  questions  of  convergence  and  error 
estimates  are  still  open.  For  example,  there  are  many  non-oscillatory  schemes  based  on 
the  minmod  limiter  which  are  numerically  robust,  at  least  in  many  numerical  tests,  but 
theoretical  results  about  convergence  and  error  estimates  are  still  missing  [3,  7,  8,  23]. 

In  this  paper,  we  consider  a  class  of  the  so-called  Godunov- type  schemes  for  solv¬ 
ing  (1).  There  are  two  main  steps  in  such  schemes:  evolution  and  projection.  In  the 
original  Godunov  scheme,  the  projection  is  onto  piecewise  constant  functions  the  cell 
averages.  In  the  general  Godunov-type  method,  the  projection  is  onto  piecewise  poly¬ 
nomials.  To  determine  the  properties  of  these  schemes  it  is  necessary  to  study  the 
properties  of  the  projection  operator.  We  limit  our  attention  to  the  case  of  piecewise 
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linear  projection  based  on  cell  averages  using  minmod  limiters  for  the  slope  reconstruc¬ 
tion  and  we  call  such  a  scheme  minmod-type.  For  example,  the  Nessyahu-Tadmor  scheme 
[19]  is  of  minmod-type  and  it  is  based  on  staggered  evolution,  other  examples  include  the 
second  order  non-oscillatory  central  schemes  with  non-staggered  grids  given  in  [9,  10]. 
Theoretical  results  about  convergence  of  such  schemes  to  the  entropy  solution,  or  error 
estimates,  are  still  missing.  In  most  cases  the  authors  give  a  variation  bound  for  such  a 
scheme  which  is  enough  to  conclude  that  the  method  converges  to  a  weak  solution.  The 
only  paper  which  has  a  convergence  result  is  the  Nessyahu-Tadmor  paper  [19]  in  which 
the  authors  prove  a  single  cell  entropy  inequality  for  a  minor  modification  of  the  original 
MinMod  scheme.  A  single  entropy  inequality  is  enough  to  conclude  that  the  scheme  is 
convergent  to  the  unique  entropy  solution  but  does  not  give  any  rate  of  convergence.  In 
order  to  get  a  rate,  one  has  to  have  a  family  of  entropy  inequalities  (see  [1,  15]).  Al¬ 
ternatively,  for  a  convex  flux,  one  can  impose  Lip+  stability  on  the  projection  and  then 
prove  convergence  via  Tadmor’s  Lip'  theory  [19,  24],  Unfortunately,  it  is  well  known  that 
minmod-type  schemes  are  incompatible  with  the  Lip+  condition  the  Lip+  seminorm 
is  not  preserved  by  a  minmod-type  projection.  It  is  easy  to  think  about  minmod-type 
schemes  in  terms  of  new/old  cell  averages.  That  is,  we  start  with  a  sequence  of  cell 
averages  {wj}  and  after  one  time  step  (projection  and  evolution)  we  get  a  new  sequence 
{w'j  }.  A  scheme  is  total  variation  diminishing  (TVD)  if  the  variation  of  the  new  sequence 
| w'j  —  w'j-i\  is  not  bigger  than  the  variation  of  the  old  one  JA  | Wj  —  wj- 1|,  i.e.,  the 
li  norm  of  the  jumps  does  not  increase  in  time.  In  the  Lip+  case,  for  convex  flux  the 
condition  on  the  jumps  is  that  the  biggest  non-negative  jump  does  not  increase  in  time 

SUp (Wj  -  (/■'_,),  <  SUp(^  -  Wj  |  ) . . 

j  j 

In  Section  3  of  this  paper,  we  prove  that  for  linear  flux  the  l2  norm  of  the  jumps  for 
some  minmod-type  schemes  does  not  increase  in  time  and  based  on  that  we  derive  a 
new  error  estimate  in  section  4.  In  the  non-linear  case,  there  is  numerical  evidence  that 
for  convex  flux  we  have  the  one-sided  analog 

Y^(wi  ~  w'j- 1)+  <  -  wj- 1)+ 

j  j 

which  suggests  a  different  approach  to  prove  convergence  for  such  schemes. 

2  Non- Oscillatory  Central  Schemes 

In  this  section,  we  are  concerned  with  non-oscillatory  central  differencing  approximations 

to  the  scalar  conservation  law 

(2)  ut  +  f(u)x  =  0. 

The  prototypes  of  all  central  schemes  are  the  staggered  form  of  the  Lax-  Friedrichs  (LxF) 
scheme  and  its  second  order  extension,  the  Nessyahu-Tadmor  (NT)  scheme  [19].  For  an 
introduction  in  central  schemes  see  [9,  10,  17,  19].  For  simplicity,  we  limit  our  attention 
to  the  staggered  NT  scheme  described  below.  Let  v(x,t)  be  an  approximate  solution 
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to  (2),  and  assume  that  the  space  mesh  Arc  and  the  time  mesh  At  are  uniform.  Let 
xj  :=  j Arc,  j  G  Z,  A  :=  At  and 

l  rxj+ 1/2 

(3)  Vj(t):=—  v(x,t)dx 

be  the  average  of  v  at  time  t  over  (rCj_i/2,  Xj+i/2).  Let  us  assume  that  v(-,  t )  is  a  piecewise 
linear  function,  and  it  is  linear  on  the  intervals  (xj^i/2,  rr,-+i/2),  j  G  Z,  of  the  form 

(4)  v(x,t)  =  Lj(x,t)  :=  Vj(t)  +  (x  —  xj)^v'j,  Xj_l/2  <  x  <  xj+l/2, 

where  fjv)  is  the  numerical  derivative  of  v  which  is  yet  to  be  determined.  Integration 
of  (2)  over  the  staggered  space-time  cell  (xj,Xj+i)  x  (t,t  +  At)  yields 


(5)  Uj+i/2(t  +  At)  = 


OC^+l/2  fxl  +  l 

Lj(x,t)dx+  \ 

Xj  j  Sj  +  1/2 


Lj+1(.x,  t)  dx 


a£-|-A£  \ 

f(v(xj+1,T))dT-  j*  f(v(xj,r))  drj  . 


The  first  two  integrals  on  the  right  of  (5)  can  be  integrated  exactly. 
CFL  condition 


(6) 


A  max  \f(y(x,  t))|  < 

Xj<X<Xj-\- 1 


1 

2’ 


j  e  z, 


Moreover,  if  the 


is  met,  then  the  last  two  integrants  on  the  right  of  (5)  are  smooth  functions  of  r. 
Hence,  they  can  be  integrated  approximately  by  the  midpoint  rule  with  third  order 
local  truncation  error.  Note  that,  in  the  case  of  zero  slopes  and  the  time 

integration  is  exact  for  any  flux  /,  and  even  for  nonzero  slopes  the  time  integration  can 
be  exact  for  a  low  degree  polynomial  flux  if  a  higher  order  quadrature  rule  is  used.  Thus, 
following  [19],  we  arrive  at 


(7)  vj+1/2(t  +  At)  -  ^(vj(t)+vj+1(t))  +  ^{v,j-v’j+1) 

~  ^(f(v(xj+i,t  +  At/2))-f(v(xj,t  +  At/2))). 


By  Taylor  expansion  and  the  conservation  law  (2),  we  obtain 

(8)  v(xj,t  +  At/2)  =  Vj(t )  -  ^A fj, 

where  stand  for  an  approximate  numerical  derivative  of  the  flux  f(v(x  =  Xj,t )). 
The  following  choices  are  widely  used  as  approximations  of  the  numerical  derivatives 
(we  drop  t  to  simplify  the  notation) 

(9)  Vj  =  m(vj+1  -vj^Vj  -Uj_i), 

fj  =  1)  -  fiVj),  f(vj)  -  f(vj- 1)), 
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where  m(a,  b )  stands  for  the  minmod  limiter 

(10)  m(a,  b)  =  MinMod(a,6)  :=  -(sgn(a)  +  sgn(fe))  •  min(|a|,  |6|) 

2 

with  the  usual  generalization 

(  inf  (E)  if  E  C  R+ 

(11)  m(E)  <  sup(F?)  if  E  C  R_  . 

^  0  otherwise 

A  generalization  of  this  numerical  approximation  is  based  the  so-called  minmod-# 
limiters 

(12)  v'  =  m  ^9(vj+ 1  -  vj),  ^(vj+ 1  -  Uj_i),  6{Vj  -  vj-^  , 

fj  =  m  (o(f(vj+ i)  -  ^(f(vj+ 1)  -  /(ni_1)),#(/(ni)  -  /(^-i))^  . 

Given  the  approximate  slopes  and  flux  derivatives  (12),  we  have  a  family  of  central 
schemes  in  the  predictor-corrector  form 

v(xj,t  +  At/2)  =  Vj(t)  -  fj, 

(13)  vj+1/2(t  +  At)  =  ~(vj(t)  +  vj+l(t))  +  g(v'j  -  v'j+1) 

~  A  (f(v(xj+1,t  +  At/2))  -  f(v(xj,t  +  At/2))) , 

where  we  start  with  Vj(0)  A-  J^3+^  uo(x)  dx.  Note  that  we  alternate  between  two 
grids  uniform  partition  of  the  real  line:  all  intervals  with  integer  end  points  for  t  =  2k At, 
k  G  Z,  and  half  integers  for  t  =  (2k  +  l)At,  k  G  Z.  As  a  special  case,  we  recover  the 
staggered  LxF  scheme  for  6  =  0  and  the  basic  MinMod  scheme  for  6  =  1  (the  middle 
slope  in  the  minmod  limiter  (12)  drops  if  6  <  1). 


3  I2  Stability  for  Linear  Flux 


In  this  section  we  will  prove  that  the  central  scheme  given  in  (13)  is  /2-stable  for  any  6  in 
the  interval  [0, 1].  Based  on  this  stability  we  will  also  derive  a  new  error  estimate  in  L2 
instead  of  the  usual  L\  estimates  in  the  conservation  laws.  Note  that  even  for  linear  flux 
/,  the  minmod- type  schemes  are  not  linear  and  the  only  global  property  known  was  that 
the  total  variation  does  not  increase  in  time  under  an  appropriate  CFL  condition,  see 
[19].  The  class  of  minmod-type  schemes  is  also  not  Lip+  stable  except  for  the  obvious 
choice  9  —  0.  Let  us  consider  a  linear  flux  f(u )  =  au,  uniform  time  steps  tn  —  nAt,  and 
restrict  the  minmod  limiter  to  9  <  1.  We  denote  :=  Vj(tn),  S'-  :=  n"  —  n”_1.  The 
minmod  scheme  (13)  reduces  to 


(14) 


?;«+i 
V3+ 1/2 


5K  +  ”5V.)  +  5 


2 
aAt 

Ax 


vn 

vj+ 1 


8 

aAt 

2Ax 


«rn(<?+1.«?+2) 


n 
j+2) 


n  ,  a At 

v  4 - 

J  2Ax 


«m  (<5”,<5"+1) 
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Hence,  we  have  an  explicit  formula  for  the  new  cell  averages  (at  time  tn+ 1)  on  staggered 
grid  in  terms  of  the  old  ones  (at  time  tn )  on  regular  grid.  In  order  to  simplify  the  notation, 
we  drop  the  time  dependence  and  denote  Wj  :=  vf,  u>'+1  :=  Sj  :=  w'-  —  , 

a  |  and  f3  |a(l  —  a).  With  this  notation,  we  have  the  following  relation 

between  the  sequence  of  the  new  averages  {w'j}  and  the  old  ones  {«</} 

(15)  w'j  —  awj- 1  +  (1  —  a)wj  +  9(3  (m(<5j_i,  Sj)  —  m (Sj,  8j+i))  ■ 

Using  that  Sj  —  Wj  —  Wj- 1,  we  derive  the  formula  for  the  sequence  of  new  jumps  in  terms 
of  the  old  ones 

(16)  8j  —  aSj- 1  +  (1  —  ot)8j  —  9(3m(8j-2 ,  <5y_i)  +  29(3m(Sj-X,Sj)  —  9(3m(Sj,Sj+%). 

The  CFL  condition  (6)  reduces  to  0  <  a  <  1  because  a  —  1/2  +  and  |^|  <  1/2. 
The  main  result  in  this  section  is  the  following  stability  result. 

Theorem  1.  If  the  initial  condition  uQ  G  L2( R),  then  the  l2  norm  of  the  jumps  of  the 
approximate  solution  is  non-increasing  in  time.  That  is 

a?)  lltPk  =  IIK+1  -  *>“-+i,}llfc  <  ll{ij}llfc  =  II K  -  ^-Jllb 

for  all  n  >  1. 

Proof.  It  is  clear  that  we  have  to  prove  the  result  for  one  step  assuming  that  |  {5/} ||/2  < 
oo.  We  split  the  proof  in  two  parts.  First,  we  prove  the  stability  for  a  monotone  sequence 
{wj}  which  is  equivalent  to  ^>0  for  all  j  G  Z.  Then  we  apply  that  result  locally  to 
derive  the  /2-stability  for  a  general  sequence. 

Theorem  2.  Let  us  assume  that  8j  >  0,  j  G  Z,  and  8j  be  given  by  (16).  Then 

(is)  ii(<pii2  <  \\m  ii2. 

Proof.  Let  us  recall  that  G  /2,  and  Sj  >  0  for  all  j.  It  is  enough  to  prove 

Theorem  2  only  for  0  <  a  <  1.  Let  (3\  9(3 ,  then  0  <  (3\  <  (3.  We  construct  the  new 
sequence  {/>'}  by  using  the  rule 

(19)  Sj  —  (1  -  a)8j  +  a8j-i  -  (3i  min (Sj-2,  Sj- 1)  +  2(3X  min(^-_i,  Sj)  -  (3X  min(Sj,  Sj+1), 

for  each  j.  First  we  assume  that  {5j}  has  finite  support.  It  is  easy  to  see  how  to  modify 
the  proof  in  case  the  support  is  not  finite.  Therefore  we  assume  Sj  —  0  for  j  <  3  and  for 
j  >  N  —  3  for  some  integer  N.  Then  Sj  —  0  for  j  <  3  and  j  >  N  —  2.  Thus  it  suffices 
to  prove 

N  N 

(20)  Z3>Z(8)2 

j= 1  j=i 

Let  us  introduce  some  notation.  Let  tjj  =  min(^-,  <5j+i),  ASj  =  Sj  —  Sj- 1,  Aijj  = 
yj  -  yj- 1,  A25j  =  Sj  -  2Sj-i  +  Sj—2 ,  and  A2y3  =  y3  -  2y3-X  +  ^_2.  Then  (19)  becomes 

Sj  =  ((1  —  a)Sj  +  OiSj- 1)  —  (3xA2yj . 
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Then 

N  N 

YM  f  =  'Y  1  _  a)^i  +  a^j~ i)2  _  2/?i((l  -  a)Sj  +  a8j-1)A2yj  +  /?i  (A2?/,-)2)  . 
j= 1  i=i 

Note  that  since  <5o  =  <5i  =  0  and  8j- ,j  =  <5j  =  0,  we  have 

N 

y~i  $j  -  ((i  - a)^ + oiSj-i)2 

3= 1 

V  V 

=  5^((!  -  (!  -  a)2)sj  -  M1  -  a)5j<5j-i  -  «2 

t=i  j=i 

N 

=  ^(1  —  (1  —  a;)2  —  a2)52  ~  2a(l  —  a)8jSj- 1 
j=:  1 

AT  N  N 

=  2,3  E^J  ~  25jdj-i)  =  2,8  Eft  -  i,,,)2  =  2/3 

3= 1  J=1 

Thus  we  get 

N  N  N 

<21)  (2/3(A<5,)2  +  2ft((l  -  «)%  +  a^-0 Aft-  -  /3?(A  V)  ■ 

j= i  i=1  i=1 

To  prove  Theorem  2,  we  need  to  prove 

N 

E  (20(A«i)2  +  2A((1  -  aft  +  <ft-i)A2»  -  f3‘^(A2yj)2)  >  0. 

3= 1 


Note  that 


T  (2/3(A<5j)2  +  2/5i((l  -  a)8j  +  a8j-l)A2yj  -  /32(A2yj)2)  — 

3= 1 

T  (2(^//?i)(A5,-)2  +  2((1  -  a)8j  +  a5j-i)A2yj  -  ^(A2^)2)^  > 

j=i 

/5i  |  ^  (2(A^-)2  +  2((1  -  a)Sj  +  a8j-1)A2yj  -  /5(A2^-)2) 

\j= i 

N 

(ft/ft  E  (2/3(A^.)2  +  2/J((l  -  a){j  +  a^Aft  -  /32(A2a)2)  . 

t= 1 

Therefore  it  is  sufficient  to  prove  the  theorem  in  the  case  f3\  =  0  (it  is  the  worst  case 
in  certain  sense). 
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Now  we  use  A 2yj  =  A yj  —  A?/j_i,  A yj  =  0,5 )  =  0  for  j  <  1,  j  >  N  —  1,  and  Abel’s 
transform  to  obtain 

N  N 

3= 1  i=1 

and 

jv  w 

j=i  j=i 

So,  (21)  becomes 

N  N 

J2S1-  E©2  =  W  (Ef=  i (A^)2  -  (1  -  a)  E"  i  Aij+1Afe 

j=i  i=i 

-a  E,E  A»  -  f  Ef=1(A2»)2)  • 

To  finish  the  proof,  it  is  sufficient  to  prove  the  following  two  Lemmas: 

Lemma  1. 

N  N  N  N 

E(A«i)2  -  (1  -  a)  E  A*i+>  A» -«E  A^A»  -  /?  E<A^)2  S  »• 

j=l  j=l  j=l  j=l 


AT  AT 

Lemma  2.  2^(A2(^)2  >  ^(A2^)2. 

j-i  i=i 

The  proof  of  Lemma  1.  We  consider  that  JT  denotes  Denote 

A  =  ^  ^  A5j.+i  Ai/j , 

j 

s  =  EAtAw. 

j 

Our  aim  is  to  prove  that 

(22)  E<At)2  -  (1  -  a) A  -  aB  -  f)  £(A%)2  >  0. 

3  j 

Let  u. |_  =  max(ii,  0),  U-  =  min(w,  0).  It  is  easy  to  check  that 

(23)  Ay,  —  (A£j)+  +  (A(5J+1)_. 

We  can  transform  A  as  follows: 

A  =  E  A<W(A<S,)+  +  (A<S,+1)-)  =  £  A^+i(Aij)+  +  EAt(At)- 

j  i  i 
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(24) 


=  ^  (A Sjf  +  J2  At  At+i 

A5j<0  ASj+1>0 

=  ^2  (^i)2  +  X/  A(5jA5j+i  +  D , 

A  5,  <0  A<5j>0,A^+i<0 

where 

D=  V  AS,A5j+1. 

A8j>0,ASj  .  i  >0 

Further, 

d  =  \  e  (w2 + (A<w2  -  <A2^+.)2) 

ASj>0,ASj+1>0 

=  \  E  «At-i)2  +  <At)2  -  <A2t)2) 

A^_i>0,A^>0 

=  \  E  <A«2+^  E  (A«2-^  E  (A^)2- 

A5j>0,A5i+i>0  A^£0,A^-_*£0  A5,-_a>0,A^>0 

By  (24), 

a  =  E(A«2  -  5  E  (A^)2-j  E  <A^)2 

i  A^>0,A^+1<0  ASj>0,ASj-i<0 

(25)  -i  E  (A2t)2+  E  AtAt«- 

A5j_i>0,A<5j>0  A^>0,A5i+i>0 

Transform  B  in  the  same  way  as  A: 

(26)  B=  E  <A^)2  +  E  'HA'5»>  +  E, 

A5j>0  A<5.,>0,A<5,,+i<0 

where 

«=  E  A^A^+t- 

A<5j<0,A5i+i<0 

The  quantity  E  can  also  be  rewritten  in  the  same  way  as  D: 

e  =  \  y  (a<5j)2  +  5  E  (A^)2  - 1  E  <A%)2- 

A5j  <0,A5j+i<0  A5j  _i<0,A5j  <0  A<5j_i<0,A^<0 

Combining  this  equality  with  (26)  we  get 

B  =  E(A«2-i  E  oh)2-^  E  (A^)2 

j  ASj<0,ASj+i>0  A8j<0,ASj-i>0 

(27)  -i  E  <A2t)2+  E  AtAt«' 

A^J-I^O.A^^O  A5,,+i<0,A(5.,>0 
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By  (25)  and  (27), 

(28)  ^(A^-)2  -  (1  -a)A-aB-p  ^(A%)2  =  F  +  G  +  H  +  I  +  J  +  K  +  L, 

j  j 

where 

p  —  —k—  E  (^)2, 

ASj  >0,  A5J+i  <0 

c  =  ~2~  E  (A«2. 

ASj>0,ASj-!<0 

H  =  %  E  (A«2. 

A^<0,A^+i>0 

i-\  E  (At)2, 

A  <  0 ,  A  Sj  |>  0 

J=-  J]  Ai,A^4.+  (i^-d  J]  (A2<5,)2 

A<5j+i<0,A<5j>0  '  '  A<5j_i>0,Ai5j>0 

+  (f  -  f>)  E  <A2^)2. 

A^_i<0,A^<0 

K  =  -/?  E  (AV. 

A5j  _|g>0,  A(53-  <0 

L  =  -/?  J]  (A2<5,)2. 

A<5j_i<0,A5j>0 

We  have  to  prove  that  F  +  G  +  H+I  +  J+K+L  >  0.  Among  the  sums  F,  G,  H,  /,  J,  K ,  L 
only  the  two  last  sums  might  be  negative;  we  will  show  that  and 

(29)  F  +  I  +  K  >  0 

(30)  G  +  H  +  L>  0. 

Indeed, 

i^A^)2  +  |(A,5,)2  -  /?(A%)2  =  ^(A^-O2  +  |(A^)2 

-  <lA5E(A^  -  A^i)2  =  i((l  -  a)A S^i  +  aAS,)f  >  0. 

Summing  the  last  inequality  over  all  j  with  ASj-i  >  0,  A5,-  <  0,  we  get  (29).  The 
inequality  (30)  can  be  proven  in  similarly. 

Also,  we  have 

(31)  J  >  0 

Finally,  plugging  (29),  (30),  and  (31)  into  (28),  we  obtain  the  required  (22).  This 
completes  the  proof  of  Lemma  1. 
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Proof  of  Lemma  2. 

First,  recall  that  A 2yj  =  0  for  j  <  1  and  j  >  N.  Also,  from  the  proof  of  Lemma  1  we 

A  8j  <  0 
A  8j  >  0 
A  5j  >  0  ' 

A  Sj  <  0 

A 5j-i  <  0 
A 5j-!  >  0 
A^_!  >  0  ' 

A<5j_i  <  0 

Therefore  A^-jj,  A  Sj,  ASJ+\  and  their  signs  determine  uniquely  A  2yj.  We  have  eight 
cases  depending  on  what  the  signs  of  A^_i,  A Sj,  ASJ+i  are. 

Case  I.  (A,  T,  T),  that  is  ASj—i  ^  0,  AS,  >  0,  A6j+1  >  0. 

Then,  A  yj  —  A  Sj,  Ayj_i  —  ASj-i,  so  A  2yj  —  A  2Sj. 

Thus,  (A 2yj)2  <  (A25j)2  in  this  case. 

Case  II.  (+,  +,  — ),  that  is  A^_i  >  0,  A  Sj  >  0,  Ah,+1  <  0. 

Then,  A  yj  —  A5J+1  +  A  Sj,  Ay^_i  =  ASj  j.  so  A  2yj  —  ASj+1  -f  A  Sj  —  ASj-i. 

Case  III.  (+,  — ,  +),  that  is  ASj-i  >  0,  A  Sj  <  0,  A5j+i  >  0. 

Then,  A  yj  —  0,  A?/j_i  =  A  Sj  +  ASj-i,  so  A  2yj  —  —  A  Sj  —  ASj^%. 

In  this  case  (A 2yj)2  —  (A 2Sj)2  =  4A<5jA<5j_i  <  0,  and  (A 2yj)2  <  (A2Sj)2. 

Case  IV.  (+,  — ,  — ),  that  is  A^_i  >  0,  A  Sj  <  0,  ASj+i  <  0. 

Then,  A  yj  =  ASj+1,  A^_i  =  A  Sj  +  ASj  i,  so  A  2yj  =  ASj+1  -  A  Sj  -  ASj-!. 

Case  V.  (— ,  +,  +),  that  is  ASj-i  <  0,  A  Sj  >  0,  A5J+1  >  0. 

Then,  A  yj  —  A  Sj,  Aijj-i  —  0,  so  A2yj  —  A  Sj. 

In  this  case  0  <  A  Sj  <  A  Sj  —  ASj-j ,  and  (A  2yj)2  <  (A%)2. 

Case  VI.  (— ,  +,  — ),  that  is  A<5j_i  <  0,  A  Sj  >  0,  ASj+i  <  0. 

Then,  A  yj  =  A<5j+i  +  A  Sj,  Ayj  ,  =  0,  so  A2yj  =  ASj+1  +  ASj. 

In  this  case  (A 2yj)2  —  (A2^-+i)2  =  4ASj+iASj  <  0,  and  (A 2yj)2  <  (A2Sj+i)2. 

Case  VII.  that  is  A5j_i  <  0,  A  Sj  <  0,  A(^+i  >  0. 

Then,  Aijj  —  0,  Ayj_i  —  A  Sj,  so  A  2yj  —  —A  Sj. 

In  this  case  0  <  —A  Sj  <  A8j+i  —  A  Sj,  and  (A  2yj)2  <  (A2<5j+i)2. 

Case  VIII.  (— ,  — ,  — ),  that  is  ASj-i  <  0,  A  Sj  <  0,  A<5,-+i  <  0. 

Then,  Ayj  =  ASj+1,  A%-_i  =  A  Sj,  so  A  2yi  =  A2Sj+1. 

In  this  case  (A 2%-)2  <  (A25j+i)2. 

Therefore  in  cases  I  (+,+,+),  III  (+,—,+),  and  V  (—,+,+),  (A 2yj)2  <  (A 2Sj)2, 
and  in  cases  VI  VII  and  VIII  (A2yj)2  <  (A%+1)2. 


AVj  = 


Similarly, 


AVj- 1 


•<SS 

i-H 

+ 

< 

ASj+1  <  0 

ASj+i  +  A  Sj  if 

ASj+1  <  0 

•<s> 

<1 

A5j+l  >  0 

0  if 

A5j+1  >  0 

(  ASj 

if  ASj  <  0 

|  ASj  +  ASj-i 

if  ASj  <  0 

Mj-i 

if  ASj  >  0 

(  o 

if  ASj  >  0 
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There  are  only  two  “bad”  cases:  II  (+,+,—)  and  IV  which  need  a  special 

treatment. 

Next,  we  define  a  sequence  of  +  and  —  signs  {s^}  where  Sj  =  +  if  A Sj  >  0  and  Sj  =  — 
if  A  Sj  <  0.  Note  that  Sj  —  +  for  j  <  3  and  j  >  N  —  2.  There  are  three  types  of  “bad” 
quadruples. 

Type  A  quadruple:  (+,+,  — ,  — ),  that  is  A^_i  >  0,  Ah,-  >  0.  A8j+i  <  0,  A<5j+2  <  0  for 
some  j. 

We  claim  that  in  this  case  the  following  inequality  holds: 

(32)  (A2Vj)2  +  (A V,)2  <  2(A 283)2  +  2(A 2Sj+1)2  +  2(A25i+2)2. 

In  this  case  A  2y3  —  Ah,-+1  +  A  Sj  —  Ah,-_i  and  A2yj+i  —  A8j+2  —  A8j  lx  —  Ah,-.  If  we 
denote  A8j-\  by  a,  A 8j  by  b,  A8j+i  by  c,  and  Ah,-+2  by  d  the  above  inequality  becomes 


(c  +  b  —  a)2  +  (d  —  c  —  b)2  <  2(6  —  a)2  +  2(c  —  b)2  +  2(d  —  c)2  for  a  >  0,  b  >  0,  c  <  0,  d  <  0, 
which  is  equivalent  to 

a 2  +  2 b2  +  2c2  +  d2  —  2 ab  +  2 ac  —  8 be  +  2 bd  —  2 cd  >  0, 
or 

(a  —  b  +  c)2  +  (b  —  c  +  d)2  —  4 be  >  0 
which  holds  since  b  >  0,  c  <  0. 

Type  B  quadruple:  (V,  ~f~,  — ,  ~i~),  that  is  A8j _ ^  h  0,A Sj  >  0,  A5J+1  <  0,  A8j+2  >  0  for 

some  j. 

We  claim  that  in  this  case  the  following  inequality  holds: 

(33)  (A2y,-)2  +  (A Vi)2  <  2(A28j)2  +  2(A2hi+1)2  +  (A2 hi+2)2- 

In  this  case  A 2yj  —  ASj+i  +  ASj  —  A6j_i  and  A2yj+i  —  A8j+2  —  Ah,-+1  —  A 8j.  Using 
the  notation  we  just  introduced  the  inequality  becomes: 


{c  +  b  —  a)2  +  (c  +  b )2  <  2(6  —  a)2  +  2 (c  —  6)2  +  (d  —  c)2  for  a  >  0,  6  >  0,  c  <  0,  d  >  0. 
Since  (d  —  c)2  >  c2  it  is  sufficient  to  prove 

(c  +  6  —  a)2  +  (c  +  6)2  <  2(6  —  a)2  +  2(c  —  6)2  +  c2,  or 

a 2  +  262  +  c2  —  2ab  +  2 ac  —  86c  >  0,  or 

(a  —  6  +  c)2  +  62  —  86c  >  0 
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which  holds  for  b  >  0,  c  <  0. 

Type  C  quadruple:  that  is  ASj  \  <  0,  ASj  >  0,  A<5,+i  <  0,  A^+2  <  0  for 

some  j. 

We  claim  that  in  this  case  the  following  inequality  holds: 

(34)  (A +  (A VO2  <  (A2Sj)2  +  2(A2Sj+l)2  +  2(A2Sj+2)2. 

In  this  case  A2yj  =  A5j+i  +  ASj  and  A2t/j+i  =  A^+2  —  ASj .  j  —  ASj.  Using  the  notation 
we  just  introduced  the  inequality  becomes: 


(c  +  b)2  +  (d  —  c  —  b )2  <  (b  —  a)2  +  2 (c  —  b)2  +  2 (d  —  c)2  for  a  <  0,  b  >  0,  c  <  0,  d  <  0. 
Since  ( b  —  a)2  >  b2  it  is  sufficient  to  prove 

(c  +  b)2  +  (d  +  c  —  b)2  <  b2  +  2(c  —  b)2  +  2 (d  —  c)2,  or 
b 2  +  2c2  +  d2  —  8 be  +  2 bd  —  2 cd  >  0,  or 
(6  —  c  +  d)2  +  c2  —  6 be  >  0 

which  holds  for  b  >  0,  c  <  0. 

We  will  call  the  type  A  (32),  type  B  (33),  and  type  C  (34)  inequalities  -  “long” 
inequalities;  we  call  the  inequalities  of  type  (A 2yj)2  <  (A2*,)2,  (A2a)2  <  (A2^+1)2  - 
“short”  inequalities. 

Now,  we  construct  a  set  of  inequalities.  We  identify  all  “bad”  quadruples  and  include 
the  corresponding  inequality  (type  A,  type  B,  or  type  C)  in  the  set.  Next,  for  all 
j  G  [1,7V]  such  that  Sj  is  not  a  middle  element  of  a  “bad”  quadruple,  and  such  that 
j  does  not  belong  to  the  “bad”  cases  II  and  IV  we  include  the  corresponding  “short” 
inequality  in  the  set.  Finally,  we  add  all  inequalities  in  the  set.  Taking  into  account 
that  A 2Sj  —  0  and  A2yJ  —  0  for  j  >  N  the  resulting  inequality  is 

N  N 

(35)  £a,(A2a)2<5>(A2,5,)2, 

j= 1  j= 1 

where  the  ajs  and  bjS  are  non-negative  integers.  To  finish  the  proof  of  the  lemma  we 
need  to  show  dj  >  1  and  bj  <  2  for  all  j  G  [1,7V]. 

Note  that  all  “long”  inequalities  have  the  form  (A 2yj)2  +  (A2yJ+i)2  <  •  •  •,  where  Sj 
and  Sj+i  are  the  middle  elements  of  a  “bad”  quadruple.  Then  a:j  >  1  if  s3  is  a  middle 
element  of  a  “bad”  quadruple.  (By  middle  element  of  a  quadruple  we  mean  second  or 
third  element  of  the  quadruple.) 

Now,  suppose  Sj  is  not  a  middle  element  of  a  “bad”  quadruple.  Then  j  does  not 
belong  to  the  “bad”  cases  II  and  IV.  Indeed  if  j  is  in  case  II:  (sj_i,  Sj ,  Sj+i)  =  (+,  +,  — ) 
then  Sj  is  a  middle  element  of  type  B  quadruple  if  Sj+2  =  +  and  a  middle  element  of  type 
A  quadruple  if  Sj+2  =  — .  Similarly,  if  j  is  in  case  IV:  (sj_i,  Sj,  Sj+i)  =  (+,  — ,  — )  then 
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Sj  is  a  middle  element  of  type  A  quadruple  if  Sj_ g  =  +  and  a  middle  element  of  type 
C  quadruple  if  .  Therefore  a  “short”  inequality  for  (A 2yj)2  has  been  included 

in  the  set  of  inequalities.  Thus  a,j  >  1  in  this  case  as  well.  We  proved  a,j  >  1  for  all 

je[i,JV]. 

Now,  we  prove  bj  <  2  for  all  j  e  [1,  AT],  Note  that  (A2^)2  can  appear  in  only  two 
“short”  inequalities:  (A 2yj)2  <  (A28j)2  and  (A2yj_i  )2  <  (A28j)2.  Therefore  bj  <  2  if 
(A2<5j)2  does  not  appear  in  any  “long”  inequalities,  that  if  Sj  is  not  a  second,  third,  or 
fourth  element  of  a  “bad”  quadruple. 

The  case  when  Sj  is  a  second,  third,  or  fourth  element  of  a  “bad”  quadruple  requires 
more  work.  First,  note  that  two  distinct  “bad”  quadruples  have  at  most  two  common 
elements.  Indeed  all  “bad”  quadruples  are  of  the  form  (*,+,—,*)  where  *  denotes  + 
or  — ,  and  no  “bad”  quadruple  has  (+,  — )  as  its  first  two  or  last  two  elements.  Next, 
the  only  case  when  two  “bad”  quadruples  have  two  common  elements  is  the  following 
configuration: 

(36)  (sj-l,  sj )  sj+ U  sj+  2,  Sj  . :{.  Sj  ,  ,1 )  =  (+,  +,—,+,  — ,  — ) 

Indeed,  Type  A  quadruple  can  not  share  exactly  two  elements  with  another  “bad” 
quadruple  because  no  “bad”  quadruple  has  (— ,— )  as  first  two  elements,  or  (+,+)  as 
last  two  elements.  Similar  analysis  shows  that  the  only  way  a  type  B  or  type  C  quadruple 
can  share  exactly  two  elements  with  another  type  B  or  type  C  quadruple  is  when  the 
configuration  (36)  occurs. 

Let  us  analyze  the  configuration  (36).  The  “long”  inequalities  which  correspond  to 
the  two  “bad”  quadruples  in  this  configuration  are: 

(37)  (A2Vj)2  +  (A2yj+1)2  <  2(A28j)2  +  2(A28j+1)2  +  (A2^)2,  and 

(38)  (A 2yi+2)2  +  (A Vs)2  <  (A2 V)2  +  2(A2 V)2  +  2(A2 5,+4)2- 

Their  sum  is 

(A2%)2  +  (A V,)2  +  (A V2)2  +  (A V3)2  < 

2(A  V2  +  2(A  V)2  +  2(A2  V)2  +  2(A2  V)2  +  2(A  2«5i+4)2- 

In  this  case  Sj,  -Sj+i,  Sj+2,  -Sj+3,  and  Sj+4  appear  as  second,  third,  or  fourth  element  of  a 
“bad”  quadruple.  Since  the  configuration  (36)  starts  with  (+,+)  and  end  with  (— ,— ) 
it  cannot  share  two  elements  with  a  “bad”  quadruple  outside  the  configuration.  This 
means  that  none  of  Sj,  Sj+i,  Sj+2,  sj+3,  and  Sj+A  can  be  a  second,  third,  or  fourth  ele¬ 
ment  of  a  “bad”  quadruple  outside  the  configuration.  In  other  words,  none  of  (A2<5j)2, 
(A2i,+1)2,  (A (A%+3f,  and  (A %+if  can  appear  in  a  “long”  inequality  other 
than  (37)  and  (38).  Since,  Sj,Sj+i,Sj+2,  and  Sj+ 3  are  middle  elements  of  “bad”  quadru¬ 
ples  (A2<5j+i)2,  (A2<V)2,  and  ( A2(iJ+3)2  can  not  appear  in  “short”  inequalities  as  well. 
Thus  bj. |_i  =  bj. |_2  =  bj+ 3  =  2.  Also,  (A 28j)2  can  not  appear  in  a  “short”  inequal¬ 
ity.  The  only  way  this  could  happen  is  (A2|/j_i)2  <  (A2<5 j)2  which  is  impossible  since 
j  —  1  is  either  in  case  I  (+,+,+)  or  case  V  (—,+,+)  depending  on  what  .Sj_2  is,  and 
in  both  cases  the  short  inequality  is  (A2?/j_i)2  <  (A2<5j_i)2.  Thus  bj  —  2.  Similarly, 
(A 28j+4)2  can  not  appear  in  a  “short”  inequality.  The  only  way  this  could  happen  is 
(A2yj+±)2  <  (A2<5j+4)2  which  is  impossible  since  j  +  4  is  either  in  case  VII  (—,—,+)  or 
case  VIII  (—,—,—)  depending  on  what  sJ+5  is,  and  in  both  cases  the  short  inequality  is 
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(A2yJ+4)2  <  (A2<5j+5)2.  Thus  bj+ 4  =  2.  This  concludes  the  analysis  of  the  configuration 
(36). 

Now,  let  Sj  be  a  second,  third  or  fourth  element  of  a  “bad”  quadruple  but  not  an  ele¬ 
ment  of  a  configuration  (36).  This  means  (A 2Sj)2  appears  in  exactly  one  “long”  inequal¬ 
ity  (it  can  not  be  a  second,  third,  or  fourth  element  of  two  distinct  “bad”  quadruples.) 
If  Sj  is  a  third  element  of  a  “bad”  quadruple  then  and  Sj  are  the  middle  elements 
of  the  quadruple  and  (A2<5j)2  does  not  appear  in  a  short  inequality.  Thus,  bj  <  2  in 
this  case.  The  cases  when  Sj  is  a  second  or  fourth  element  of  a  “bad”  quadruple  need 
separate  consideration. 

1)  Sj  is  a  second  element  of  type  A  quadruple  (+,  +,  — ,  — ).  The  only  way  (A2<5j)2 
could  appear  in  a  “short”  inequality  is  (A2?/j_i  )2  <  (A 25j)2  which  is  impossible  since 
j  —  1  is  either  in  case  I  (+,  +,  +)  or  case  V  (— ,  +,  +)  depending  on  what  Sj_2  is,  and  in 
both  cases  the  short  inequality  is  (A2?/j_i)2  <  (A2^_i)2.  Thus  bj  =  2. 

2)  Sj  is  a  fourth  element  of  type  A  quadruple  (+,+,—,—).  The  only  way  (A25j)2 
could  appear  in  a  “short”  inequality  is  (A 2yj)2  <  (A2Sj)2  which  is  impossible  since  j  is 
either  in  case  VII  (—,—,+)  or  case  VIII  (—,—,—)  depending  on  what  Sj+ 1  is,  and  in 
both  cases  the  short  inequality  is  (A 2yj)2  <  (A28J+\  )2.  Thus  bj  =  2. 

3)  Sj  is  a  second  element  of  type  B  quadruple  (+,+,—,+).  Here  the  argument  is 
word  by  word  like  in  1).  The  only  way  (A2Sj)2  could  appear  in  a  “short”  inequality 
is  (A2?/j_i)2  <  (A 2Sj)2  which  is  impossible  since  j  —  1  is  either  in  case  I  (+,+,+)  or 
case  V  (—,+,+)  depending  on  what  Sj- 2  is,  and  in  both  cases  the  short  inequality  is 
(A VO2  <  (A28j  ,)2.  Thus  bj  -  2. 

4)  Sj  is  a  fourth  element  of  type  B  quadruple  (+,+,—,+).  Since  the  coefficient  of 
(A2<5 j)2  in  the  corresponding  “long”  inequality  (33)  is  1  and  (A2<5j)2  could  appear  in  at 
most  1  “short”  inequality  bj  <  2. 

5)  Sj  is  a  second  element  of  type  C  quadruple  (— ,  +,  — ,  — ).  Since  the  coefficient  of 
(A28j)2  in  the  corresponding  “long”  inequality  (34)  is  1  and  (A2SJ)2  could  appear  in  at 
most  1  “short”  inequality  bj  <  2. 

6)  Sj  is  a  fourth  element  of  type  C  quadruple  (—,+,—,—).  Here  the  argument  is 
word  by  word  like  in  2).  The  only  way  (A2^)2  could  appear  in  a  “short”  inequality 
is  (A 2yj)2  <  (A 28j)2  which  is  impossible  since  j  is  either  in  case  VII  (—,—,+)  or  case 
VIII  (  —  ,—,—)  depending  on  what  s,+i  is,  and  in  both  cases  the  short  inequality  is 
(4W  <  (AV)».  Thus  bj  =  2. 

We  have  shown  that  in  all  cases  bj  <  2  for  j  G  [1,  N]  which  completes  the  proof  of 
Theorem  2.  □ 

Now,  we  continue  with  the  general  case.  That  is,  we  want  to  show  that  the  Z2  norms 
inequality 

(39)  IIU'}||i2<ll{^}||«2 

holds  for  any  initial  sequence  {^  }  with  finite  /2  norm.  We  consider  the  sequence  {wj} 
and  restrict  the  index  j  to  a  maximal  subset  Am  on  which  the  piecewise  constant  function 
w  is  monotone,  recall  that  Sj  —  Wj  —  Wj-\.  Given  a  sequence  {wj},  we  can  decompose 
it  into  monotone  subsequences.  This  decomposition  also  gives  a  decomposition  of  the 
sequence  {<5j}  into  nonnegative  subsequences.  Without  any  limitations,  we  assume  that 
the  jumps  {^  }  are  non- negative  for  all  l  <  j  <  r*  Si- 1  <  0  and  <5r+i  <  0.  That  is,  1 
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ia  a  local  minimum  and  wr+1  is  a  local  maximum  of  the  piecewise  constant  function  w. 
Let  wm  be  the  following  piecewise  constant  correction  of  w 

{Wj,  if  l  <  j  <  r, 

wi- 1,  if  j  <  l , 

wr+l,  if  j  >  r. 

<  r  +  1}  and  the  jumps  sequence  6m  {8™}  of  wm  is  given 

_  /  wj  ~  wj~  i)  if  /  <  j  <  r  +  1, 
f  0,  otherwise 

Hence,  we  have  a  sequence  of  monotone  functions  {wm}  and  the  corresponding  jump 
sequences  such  that 

E  E  iift”}ii>2  =  E  iK^niife  =  iift}ii.2 

m  jeAm  m,j£Z 

because  the  sequence  of  the  jumps  of  {ci,}  is  decomposed  into  disjoint  jump  subsequences 
{(jj"}.  There  are  two  types  of  jumps  <5'.  A  jump  Sj  is  of  type  1  if  it  is  equal  to  the  jump 
that  is  the  jump  generates  with  the  starting  sequence  {<5"'},  where  the  index  m 
such  that  j  G  Am.  A  jump  is  of  type  2  if  it  is  not  of  type  1.  Note  that  a  type  2  jump 
5j»  occur  only  inside  an  interval  which  contains  a  strict  local  extremum.  Near  a  local 
extremum  we  have  two  nonzero  jumps,  say  5ljt  and  5'-, ,  generated  by  the  two  monotone 
wm-s  with  index  sets  finishing/starting  with  j*.  It  is  easy  to  verify  that 

K-IHI4I  —  M- 

Hence,  we  have  that 

(^)2<(^)2  +  (^)2, 

and  we  conclude  that 


(40) 

Note  that  Am  —  {j  :  l  <  j 
by 

(41)  5f 


<3)2  <EE  <3(^))2  <  E  E  <4)2  =  Eft)2. 

j  m  jeAm  m  je Am  n 

where  we  use  the  notation  8j(5m)  for  the  new  jumps  generated  by  8rn.  It  is  also  easy  to 
prove  a  local  inequality  but  with  index  set  for  <5'  starting  from  an  interval  right  after  an 
extremum  and  finishing  right  before  one.  □ 


4  Error  Estimates  for  Linear  Flux 

Recall  that  u  is  the  entropy  solution  to  the  conservation  law  ut  +  f(u)x  —  0  with  initial 
condition  u°,  and  v  is  the  the  numerical  solution  described  in  (13).  In  the  case  of  linear 
flux  and  0  <  9  <  1,  the  formula  for  the  new  averages  of  the  minmod  scheme  is  given  in 
(14)  and  the  conservation  laws  (2)  reduces  to 

(42)  ut  +  aux  =  0. 


15 


Let  ST  be  the  shift  operator  defined  by  STg(-)  :=  g( ■  —  r).  Then  the  exact  solution 
of  (42)  at  time  t  for  any  initial  data  u°  is  u(t ,  •)  =  Satu°.  Let  Ah  be  the  averaging 
operator  defined  on  a  uniform  partition  by  Ahg\i  :=  ^  Jf  g(s)ds,  where  |/|  =  h.  It  will 
be  useful  to  define  a  global  approximate  solution  v.  We  first  define  the  approximate 
solution  at  discrete  times  by  vn  v(nAt,  •),  n  —  0, 1 . . . ,  N,  in  the  following  way:  (i) 
v°  :=  u°;  (ii)  vn  :=  SaAtPhvn~1 ,  where  for  odd  n,  1  <  n  <  N,  PhV  is  the  linear  function 
on  Ij  :=  (xj_1/2,Xj+i/2)  defined  in  (4)  with  the  minmod  slopes  (9),  and  for  even  n  we 
have  the  analog  definition  of  Ph  on  the  shifted  partition  {/j+i/2|  j  G  Z}.  Note  that, 
Phvn  —  PhAhVn  because  the  piecewise  linear  projection  P/,,  defined  in  (4)  and  (9),  is 
based  only  on  the  averages  of  vn  on  the  corresponding  partition.  The  formula  (14)  for 
the  new  cell  averages  can  be  written  as 

””+1/2  =  Mv")  l/J+1/1  =  MSimPkv"^)  |,j+1„ 

for  odd  n,  with  Ah  based  on  the  staggered  partition  {/j+i/2|  j  G  Z}  and  Ph  based  on 
regular  partition  {Ij\j  G  Z}.  For  even  n,  we  have  the  same  sequence  of  operators  but  on 
the  reversed  partitions.  The  global  approximate  solution  v  is  defined  by  v(nAt,  •)  =  vn 
and  v(t,  •)  =  Sa(t_n/±t)(Ph,vn)  for  nAt  <  t  <  (n  +  1) At  and  n  —  0, 1, . . . ,  N  —  1.  That  is  v 
solves  exactly  (42)  for  nAt  <  t  <  {n  +  l)Af  with  initial  data  PhVn,  n  =  0, 1, . . . ,  N  —  1. 

In  order  to  describe  the  next  result,  we  need  to  introduce  some  notation.  A  function 
g  is  of  bounded  variation,  i.e.,  g  G  BV(R),  if 

n 

I g | b v (r)  :=  swp^2\g(xi+1)  -  g(xi)\  <  oo, 

i=i 

where  the  supremum  is  taken  over  all  finite  sequences  X\  <  . . .  <  xn  in  R.  Functions  of 
bounded  variation  have  at  most  countable  many  discontinuities,  and  their  left  and  right 
limits  g{x~)  and  g(x+)  exist  at  each  point  x  G  R.  Since  the  values  of  the  initial  condition 
u°  on  a  set  of  measure  zero  have  no  influence  on  the  numerical  solution  v  and  the  entropy 
solution  solution  u,  it  is  desirable  to  replace  the  seminorm  |  •  |bv(r)  by  a  similar  quantity 
independent  of  the  function  values  on  sets  of  measure  zero.  The  standard  approach  in 
conservation  laws  is  to  consider  the  space  Lip(l,L1(R))  of  all  functions  g  G  L1(R)  such 
that  the  seminorm 

(43)  MLiP(i,Li(R))  :=  bm sup  -  /  | g(x  +  s)  -  g(x)\dx 

is  finite.  It  is  clear  that  |5f|LiP(i,L1(R))  will  not  change  if  g  is  modified  on  a  set  of  measure 
zero.  At  the  same  time  the  above  two  seminorms  are  equal  for  functions  g  G  BV(R) 
such  that  the  value  of  g  at  a  point  of  discontinuity  lies  between  g(x~)  and  g(x+)  (see 
Theorem  9.3  in  [6]).  Similarly,  we  define  the  space  Lip(l,  LP(R)),  1  <  p  <  oo,  which  is 
the  set  of  all  functions  g  G  Lp( R)  for  which 

(44)  ||#(-  -  s)  -  p(-)||lp(r)  <  Mh ,  5  >  0. 

The  smallest  M  >  0  for  which  (44)  holds  is  |g|Lip(i,LP(R))-  It  is  easy  to  see  that  in  the 
case  p  —  1  the  seminorm  given  in  (44)  is  the  same  as  the  one  in  (43).  In  the  case  p  >  1, 
the  space  Lip(l,Lp(R))  is  essentially  the  same  as  VF1(LP(R)),  see  [6]  for  details.  With 
this  notation,  we  have  the  following  result. 
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Theorem  3.  Let  u  be  the  entropy  solution  to  (2)  with  linear  flux  f(z )  =  az  and  v  be 
the  numerical  solution  described  in  (14)  with  0  <  6  <  1.  If  the  CFL  condition  (6)  is 
satisfied,  tn  =  nAt,  0  <  n  <  N,  and  T  =  N At,  we  have 

(45)  | \u(T,  ■)  -  v(T,  OII^r)  <  C'(A/’/i)1'/2/i1/2|'U°|Lip(1)LP(K)), 
for  p  =1,2  where  C  is  an  absolute  constant. 

Proof.  The  L\  estimate  is  based  on  the  TVD  property  of  the  numerical  solution  v  and 
the  L-2  estimate  is  based  on  the  l2  stability  of  the  jumps  proved  in  Theorem  1.  Both 
estimates  use  a  dual  argument  similar  to  the  one  in  [24]  and  in  the  proof  we  use  an  index 
p ,  where  p  G  {1,2}.  In  the  proof,  C  will  be  an  absolute  constant  that  can  be  different 
at  different  places. 

Let  eft,x)  :=  u(t,x )  —  v(t,x),  and  Eft,  x)  :=  ff  e(t,s)ds ,  where  we  assume  that 
u°  G  L1(W)  to  guarantee  that  E  is  well  defined  for  all  ft,  x)  G  (0,  T)  x  R.  We  have  that  E 
also  satisfies  (42)  for  nAt  <  t  <  (n  +  l)Af  with  initial  data  ff  ufnAt,  s)  —  PhVn(s)ds, 
n  =  0,1 , ...  ,N  —  1.  For  a  function  g  G  L1(R),  we  define  a  minus  one  norm  in  the 
following  way 

(46)  :=  II  /  g(s)ds\\LP{R). 

J  —  oo 

It  is  easy  to  verify  that  for  any  r  G  R 

(47)  HXsIl-^p  =  ||ff  11-i.p- 

Recall  that  T  =  NAT.  Then,  we  have  the  representations  ufT,  •)  =  ( SaAt)Nu°  and 
v(T,  •)  =  ( SaAtPh)Nu° ■  Using  (47),  we  have 

II e(T,  -)||-i,P  =  \\(SaAt)Nu°  -  (SaAtPh)Nu° ||_1|P  =  IKSU^-V  -  PhiSaAtPh)"-1*0 
and  by  the  triangle  inequality  we  obtain 

(48)  ||e(r,.)||-i„  <  ||(5aAt)JV-1u0-(^AtP,)JV-1u0||_1;P 

+  ||n(^A/)  V”,«0  -  (PaAtPh)JV“1«°||-i,P. 

Let  en  —  ((SaAt)n  ~  ( SaAtPh)n )  u°,  n  —  0, 1, . . . ,  N.  Then  (48)  is  equivalent  to 

(49)  ||eJV||_llP  <  ||eJV_1||_1)P  +  H^"1  -  nJV-1||_1)P, 
and  applying  (49)  for  n  —  N,  N  —  1, . . . ,  1,  we  get 

N 

(50)  ||eJV||_1;P<^||P^-u”||_1,p 

n=l 

because  e°  =  0.  To  prove  the  error  estimates,  we  need  the  following  technical  lemma. 
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Lemma  4.  For  any  p  G  {1,  2}  and  any  n  =  0, 1, . . . ,  N,  we  have 


(i)  IK^jllk  <  h1  p||«°||Lip(l,LP(R)), 

(ii)  (yn)’ ftl+’lltt}lk,. 

(iii)  \\Ahvn  <  f^fj)  A2|l“°llLip(i,LP(ii))- 


Proof.  The  inequalities  (i)  and  (ii)  follow  by  standard  arguments,  therefore  we  only  prove 
(i)  in  the  case  p  =  2  and  omit  the  rest.  Recall  that  8f  =  vj  —  and  by  Theorem  1 
we  have 


5 


where  8®  =  uQ-  —  u°-_ l,  u®  :=  ^  ff  u°(s)ds.  Hence,  to  prove  (i)  for  p  =  2,  we  need  to  prove 

<  ^ll'u°lll,ip(l,L2(R))- 


Since 


\s  +  h)  —  u°(s))<isj  <  h  2  yj  lM°(s  +  h)  —  n°(s)|ds  J  , 
and  since  by  Cauchy-Schwartz 

\  2 

|tt°(s  +  h)  —  u°(s)|d,sj  <hj  |u°(s  +  h)  —  u°(s)\2ds, 

we  obtain 

(51)  ^(<S°)2  <  h_1  f  |tt°(s  +  h)  —  n°(s)|2ds. 

j 

^From  (44),  we  have 

f  \n\s  +  h)-u\s)\2ds<h2\\nYup(iW)) 

J  R 

and  using  that  in  (51),  we  conclude 


—  ^■||'u°llLip(l,L2(R)) 

j 

which  proves  (i)  for  p  —  2.  To  prove  (iii),  we  note  that 

(52)  \AhVn  —  vn\ /.  <  m.SLXVn(x)  —  minnn(rr), 

3  X^Ij  X^Ij 
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and  because  vn  =  SaAt;un  1  we  have  that 

max/(a;)  —  minn"(:r)  <  2 max(|<^Ti  |,  |5”_1|). 

X(^Ij  X  £  1 j 

The  rest  of  the  proof  of  (iii)  is  analogous  to  the  proof  of  (i).  □ 

Combining  (i)-(iii),  we  have 

\\PhVn  —  nn||_iiP  <  Ch2||u0||Lip(i)Lp(R)) 

and  after  applying  the  above  inequality  in  (50),  we  derive  the  following  estimate 

(53)  lleiV II  — l,p  <  C Nh? ||u° ||Lip(l;LP(K))  • 

Because  vN  ^  Lip(l,  L2(M)),  we  approximate  vN  by 

I  nx+h/2 

V  :=  -  AhvN(s )  ds. 

II  Jx-h/2 

Similar  to  Lemma  4,  it  is  easy  to  verify  that  for  p  G  {1,2}  we  have 

(54)  P  -  vN\\-i,p  <  C'h2||w0||Lip(i)Lp(K)), 

(55)  ll^  —  WJV||lp(R)  <  C'h||w°||Lip(l,LP(R))) 
and 

(56)  P||Lip(l,LP(R))  <  P°||Lip(l,LP(M))  • 

Let  e  :=  u(T,  •)  —  v.  Then 

II  ~ll  ^  ||  JV 1 1  ,  ||  ~  JV  n 

II e  1 1  — i,p  —  1 1 e  ||-i,p+P  —  v  ||-i,p) 
and  combining  the  estimates  (53)  and  (54),  we  have 

(57)  Pll-yp  <  C  Nh?  1 1  "li0 1 1  Lip(l,LP  (®))  • 

Kolmogorov-Landau  inequalities  in  LP(W)  (page  156  in  [6])  for  the  functions  E(x)  :  — 
foo  e(s)  ds,  E' ,  and  E"  give 

IIB'IU,  <  V2||B||t,2||B"||i/»2, 

and  using  (57)  and  (56)  we  arrive  at 

(58)  ||e||iP  <  C'(A^h)1//2h1/2||w0||Lip(i,LP(K))- 
Finally,  by  the  triangle  inequality 

||eIUp  <  ||e|| lp  +  ||#  _  wJV||lp(r)  <  C,(A/’h)1'/2h1'/2||u0||Lip(i)Lp(K)), 

and  we  combine  (55)  and  (58)  to  conclude 

||«(r,.)-«(T,.)||  LP(M)  —  1 1 ^ 1 1 LP  <  C(iV/l)1/2/l1/2||tt«||Lip(1;LP(R)). 

Note  that  C  can  be  computed  explicitly  and  it  is  not  very  big  (C  <  20).  In  the  case 
p  =  2  and  u°  L1(R),  we  get  the  same  error  estimate  via  an  approximation  procedure 
because  the  estimate  is  independent  of  the  L 1  norm.  □ 
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Corollary  5.  In  the  case  of  Nh  <  C ,  we  get  the  convergence  rate 


U(T,-)  —  v(T,  -)HiP(K)  <  C'/?.1//2|u0|Lip(l,LP(K)), 


for  p  —  1  and  p  —  2. 

The  L 1  estimate  is  not  new,  it  follows  from  the  arguments  in  [24],  but  the  1/2  rate 
in  L2  is  new.  Note  that,  using  the  L1  estimate,  by  interpolation  arguments,  we  get  only 
1/4  rate  in  L2. 
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